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Translating the dynamics of the Henon-Heiles (HH) hamil- 
tonian as a geodesic flow on a Finsler manifold, we obtain a 
local and synthetic Geometric Indicator of Chaos (GIC) for 
two degrees of freedom (dof) continuous dynamical systems 
(DS's). It represents a link between local quantities and 
asymptotic behaviour of orbits and gives a strikingly evident, 
one-to-one, correspondence between geometry and instabil- 
ity. Going beyond the results attained using the customary 
dynamical approach and improving also on the global crite- 
ria established within the Riemannian framework, the GIC is 
able to discriminate between regular and stochastic orbits on 
a given energy surface, simply on the basis of the value it as- 
sumes along a relatively small piece of the trajectory, without 
long integrations of the dynamics and without any reference 
to a perturbation. 

PACS: 05.45. +b ; 02.40.-k ; 02.90.+p 



The presence of instability in DS's is usually recog- 
nized a posteriori, looking at the long time evolution of 
small perturbations, whose average exponential growth 
is interpreted as the signature of Chaos. Skipping here 
most of the issues related to the universal meaning at- 
tributed to it, we remark how the procedures and tools 
having at the grounds such a criterion alone manifest 
their limitations whenever is investigated the behaviour 
of DS's at the boundary between quasi-integrability and 
stochasticity (e.g., Jl|,|^|). Instead, the search for a priori, 
synthetic signatures of Chaos, dates back to Toda || and 
continued, across interesting investigations on the mech- 
anisms of transition to stochasticity (see e.g. Q), up 
to the recently revived Riemannian Geometrodynamical 
approach (GDA) [^J. Though most of its results relate 
to high dimensional hamiltonian systems (for which some 
approximations are justified by the large number of dof or 
some weak form of the ergodic hypotesis), more recently, 
this approach has been tested also for small DS's, giving 
outcomes clearly supporting its reliability, Never- 
theless, if in the case of large DS's considered the agree- 
ment between the GDA and the customary tools used to 
detect Chaos has revealed to be rather satisfactory, at 
least as long as the approximations are well justified 
some discrepancies emerge in the case of few dof systems 
For the latter, the GDA provided an alternative 
way to recover most of the results obtained with the tan- 
gent dynamics equations, suggesting deeper hints for the 
understanding of the sources of Chaos and giving in ad- 



dition some global criteria to single out a transition in 
the overall behaviour 0. However, this criterion is un- 
able to correctly detect the occurrence of Chaos in single 
orbits Jio| , as it renounces, in principle, to intrinsically 
describe the behaviour of individual trajectories. Within 
the Riemannian approach to few dof systems, this is- 
sue has been addressed, up to now, only resorting to 
a numerical procedure analogous to the integration of 
the tangent dynamics equations, whose results have been 
generally confirmed (though not always). Recently Kan- 
drup JnJ investigated in details the relationships existing 
between local dynamical behaviour and local geometric 
features of the Jacobi (Riemannian) manifold for some 
two-dimensional DS's, obtaining qualitative correlations 
among average curvature and its fluctuations and some- 
what more ambiguous ones between curvature fluctua- 
tions and short time Lyapunov exponents. In summary, 
the Riemannian GDA has been able, up to now, either to 
intrinsically describe the average behaviour of a DS or to 
single out the individual orbits instability a posteriori, as 
in the Hamiltonian description. Both these approaches 
have revealed unable to find an intrinsic and individual 
indicator of long-term behaviour of orbits as the GIC here 
presented, built within the Finsler GDA and which cannot 
even to be defined, (for the H H system) within the other 
above mentioned frameworks. Notwithstanding its local 
character (in both spatial and temporal meanings), it re- 
veals to be unambiguously related to Lyapunov Charac- 
teristic Numbers (LCN's), i.e., to asymptotic quantities, 
usually computed with reference to a perturbation! We 
claim then that it represents a strong indication (if not 
a proof) that the GDA is able not only to reproduce and 
to explain the results obtained with the usual tools, but 
even to go beyond them. 

One of the main recent results of the GDA is the confir- 
mation that the onset of unpredictability in the geomet- 
ric transcription of realistic (large) DS's is driven by the 
mechanism of parametric instability ]5|,^| and thus dif- 
fers completely from what occurs in the geodesic flows of 
abstract Ergodic theory. Indeed, most of phase space of 
large physical DS's is not characterized by (constant) neg- 
ative curvatures, but stochasticity is caused by the quasi 
random fluctuations of (mostly) positive curvatures. For 
few dimensional DS's however, such random character 
cannot be assumed and instead is parametric resonance, 
similar to that occurring in the Mathieu equation, to 
bring about instability, Pj7p^|. Nevertheless, we find 
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below that even for two dof systems such a mechanism 
cannot be singled out in a naive way and very intricate 
combinations of geometric features of the manifold are 
linked to instability in a non trivial way. 
We already discussed the motivations for an extension 
of the GDA to include non Riemannian manifolds JI^ ] 
and we also pointed out its somewhat greater effective- 
ness with respect to the usual tools in the computations 
of instability exponents |2|,|l3|]. We refer to ||||D f° r a 
detailed description of the GDA, to P, p"2|jL4| for its imple- 
mentation within Finsler manifolds and to [p| for a thor- 
ough discussion of most of the points here only sketched. 
Within the GDA the trajectories of a AA dof system 
become the geodesies of suitable differential manifolds, 
which, in Finsler geometry, are (TV + 1)— dimensional and 
represent a generalization of Riemannian ones. The sta- 
bility of the flow is determined by the Jacobi-Levi-Civita 
(JLC) equations of geodesic spread: 
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ds 
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H a c z c = , (a = 0,1,..., AT) , (1) 



being z a the perturbation, V/ds the covariant deriva- 
tive along the geodesic and the stability tensor, Ti., H|), 
derives from the (generalized) curvature tensor of the 
manifold. The Finsler time-parameter s = sf is defined 
through the Lagrangian function C, as dsF = Cdt and 
possesses a built-in invariance with respect to an arbi- 
trary rescaling of the Newtonian time t, jl^]. The local 
behaviour of geodesies is determined by the eigenvalues 
of Ti, which are the principal sectional curvatures (psc's) 
defined by the given geodesic on the manifold ||. For 
a A/" degrees of freedom DS, once a geodesic is chosen, 
the Finsler stability tensor possesses (TV + 1) eigenval- 
ues, ({Aa}, A = 0, 1, . . . , A/") one of which vanishes iden- 
tically, Aq = 0, associated with a neutral eigenvector, 
along the tangent to the geodesic. In the case of a stan- 
dard hamiltonian system (i.e. without gyroscopic terms), 
H{q_,p) = | a l jp t p j +U(q) we have then, 



Aj = t'(B + t' pi 



(i 



1. 
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(2) 



where B is related to the time derivatives of the La- 
grangian, the {pi} are the eigenvalues of the hessian 
U,ij (q) and the prime denotes the derivative wrt to sp- 
in this letter we will deal with the well known two di- 
mensional HH system, whose Hamiltonian is 

H( q ,p) = ±(x 2 + y 2 +pl+p 2 y )+x 2 y-^. (3) 

We find that for this DS (as well as for generic realistic 
ones, in spite of some persistent claims, e.g. |l(|) nega- 
tive curvatures are quite unable to explain the asymptotic 
character of orbits. The mechanism of parametric insta- 
bility, due to fluctuations of usually positive curvatures 
^,[1 manifests naively however only for some many dof 



systems, being instead hardly perceived in this case. For 
example, the analysis of spectra, ||], shows how inter- 
mingled and far from trivial are the relationships with 
the elementary theory of Mathieu-like equations. Such 
seemingly discouraging inconsistencies driven us to check 
for different signatures of instability. The Finslerian ap- 
proach allows to consider, even for two dof Hamiltonians, 
the possible anisotropy of the manifold, which results to 
play a crucial role, via the Schur theorem Q in the mech- 
anism of instability: fluctuating curvatures require also 
that the manifold is anisotropic. The connection between 
curvatures variations along a geodesic and anisotropy, on 
one side, and growth and rotation of perturbation, on 
the other, is far from being clear and is currently investi- 
gated: how the former can interact to steer the latter can 
at the moment only be guessed. For a two dof system, the 
associated Finsler manifold is three dimensional and its 
curvature properties along a given geodesic are described 
by the two non vanishing psc's, Ai j2 , which are invariant 
functions on the tangent space, representing the sectional 
(i.e. gaussian) curvatures in the two-planes defined from 
the flow and the two (non trivial) eigenvectors of Ti. 
Given them, we can characterize the way the geodesic 
explores the manifold through the (half) Ricci curva- 
ture (along the flow) and the anisotropy, «[q(s),p(s)] 
and t9[q(s), p(s)], respectively: 

. Ai + A 2 _ Tr(W) Ric F (u) „ . Ai - A 2 







(4) 



An exhaustive statistical analysis of the behaviour of k 
and -d (or equivalently of Ai, 2 ) and the details of the log- 
ical path leading to the synthetic indicator are presented 
elsewhere, ||. We found that along a generic geodesic 
the {Xi} oscillate around their average values, the fluctu- 
ations of Ricci curvature in general turns out to be how- 
ever much smaller than those of the sectional ones, which 
are indeed almost anti-correlated. Such an effect is par- 
ticularly evident in the HH ISLA — 2. So, the man- 
ifold appear to be everywhere anisotropic but with psc's 
always (Ai) or mostly (A 2 ) positive. Fluctuations (and 
then anisotropy) increase with energy, showing a global 
qualitative change in correspondence of appearance of 
stochasticity, [||. Although smaller, the overall fluctu- 
ations of Ricci curvature seem nevertheless to influence 
appreciably the stability of geodesies. Moreover, Schur 
theorem asserts that the two quantities must be related. 
We then look for a relative measure of anisotropy fluctu- 
ations compared to overall curvature variations. Corre- 
lations between two quantities A(s) and B(s) reflect in 
the functional: 

C S [A,B] = {A ' B)s (5) 

((A 2 )s ■ (B 2 ) s f 2 

which clearly depends on the averaging interval S, depen- 
dence understood in the sequel. An indication about the 
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relative importance of anisotropy wrt average curvature 
fluctuations can be obtained comparing the phase space 
normalized correlation functions C[i9, 5$] and C[k,5k], be- 
ing SA(s) = A(s) — (A). Using them we build up a quan- 
tity which apparently describes only the local geometric 
features of the sub-manifold explored (in a lapse S) by the 
geodesic having initial conditions [q(0),p(0)] = (q ,p ): 

R F [S] = i?F(q ,P ) = §M ^ ° ^ 

C[K, 0K\ 

However, the inspection of figure [l] shows a striking corre- 
spondence between the PSS's and the map of R 2 F . Indeed, 
figures H a) and b) , which refer to the typical energy value 
E = 1/8, show that smaller is the value of Rp(q ,p ), 
more regular is the geodesic passing through (q , p ). To 
obtain the plot of figure [j] b) , we chosen a grid of points 
{y ,p Vo } on the PSS x — 0, chosingp Xo such that H = E 
and then numerically integrated the geodesic equations, 
computing the correlation functions entering i?F(q , P )- 
The same results have been obtained at all the energies; 
the more essential, though less suggestive, histograms of 
figure ||, confirm that our GCI is able to depict correctly 
the single orbits behaviour up to the dissociation energy. 
The columns in this plot represent the values of R 2 Fe , 
defined as 

fl| E (q ,Po) = a 2 (E) #!(q ,Po) , (7) 

where a(E) — a± ■ E 2 is a scaling factor to get rids of 
a global energy dependence, useful to compare results at 
different energies. For each energy, the R 2 F values are 
reported for a sample of seven initial conditions, chosen 
among those topologically equivalent to the ones indi- 
cated in figure |l|b) for E = 1/8. We observe that while 
chaotic trajectories are always characterized by Rp„ val- 
ues around its upper limit (normalized to unity this is 
because we use R 2 p E instead of Rf), regular orbits have 
instead considerably smaller values, which is however 
higher for those geodesies tending to become chaotic ear- 
lier, as the energy increases. In particular, from both the 
figures, we see that the geodesies in the regular islands 
located on the p y — axis of the PSS are recognized by our 
GIC as nearer to chaotic ones than those belonging either 
to the large regular island on the y— axis or to the banana 
region; this explains why these islands disappear earlier 
as the energy increases and gives also some insights on the 
causes of the (partial) failure of the Riemannian approach 
to describe these peculiar orbits P,|7|]Tc|]. Moreover, very 
interesting insights can be obtained looking at the relax- 
ation patterns of Rf as function of the averaging time 
S P], as can be perceived from the diffusive behaviour 
around the border of regular islands. Well then, a ge- 
ometric local quantity turns out to be deeply related to 
the asymptotic behaviour of geodesies. We stress that the 
values plotted in the map and histograms are obtained 
through computations of correlation functions over time 



intervals much shorter than those needed to obtain either 
the PSS's or the LCN's: it suffices to follow a geodesic for 
as few as fifteen periods in order to see if the value of 
Rf attains the upper limit characterising the stochastic 
sea or tends to a lower value, which though different for 
distinct regular orbits, is however always smaller than for 
chaotic ones. A local signature of asymptotic instability 
acquires a special significance also on the light of the is- 
sue of reliability of long time numerical integrations of 
chaotic systems [jl8). The definition of Rp obviously im- 
plies that it can be used also as a global indicator, able 
to give a quantitative measure of the overall degree of 
stochasticity at a fixed value of energy: a phase-averaged 
local indicator is a global one, instead a global indicator 
is, in general, non-local. Perplexities which can be raised 
by the apparently cumbersome definition of Rp are an- 
swered on the light of the pathologies affecting the HH 
hamiltonian, which amount mainly to the degeneracies 
of its integrable limit (TJJ, p. 46]. Indeed, for most two 
dimensional DS's, more naive Finsler geometric indica- 
tors suffice to discriminate between chaotic and regular 
orbits. Although the outcomes here presented (and in 
clearly support the reliability of Finsler GDA, nev- 
ertheless a better understanding is still waiting, as we 
need to test and extend the proposed criterion to more 
general and higher dimensional DS's. Moreover it is de- 
sirable to improve the rather phenomenological interpre- 
tation of Rf and possibly to predict theoretically the 
threshold above which stochasticity occurs. This goal 
amounts essentially to understand whether high Rp val- 
ues either follow from or are a cause of stochastic be- 
haviour (or both). The results obtained support the con- 
viction mm that negative curvature is unnecessary (if 

not irrelevant) to explain the onset of Chaos in realistic 
DS's, keeping however in mind that some faded global 
correlations exist. Instead, a separate spot concerns the 
long-lasting claims about the implications of scalar curva- 
ture: except for two dimensional or isotropic manifolds, 
it appears to have nothing to do with the behaviour of 
geodesies representing realistic DS's || and an hopefully 
coherent explanation of its irrelevance will be presented 
elsewhere. Among the issues still open, it is worth to 
investigate more deeply the relevance of the rotation of 
eigendirections oiH: a direct inspection of the JLC equa- 
tions seems to indicate that the rotation of the pertur- 
bation vector should give a contribution to instability, 
nevertheless a separation by hand of interacting effects 
often causes serious inconsistencies. 
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FIG. 1. Phase space portraits of HH system at E = 1/8: 
a) PSS obtained numerically integrating a sample of orbits 
up to T— 10000 units; b) Gray-scale plot of RF(y ,p yo ) com- 
puted following a set of geodesies for a lapse 5=200 units. 
Darker dots represent smaller Rf values, the White Stochas- 
tic Sea corresponds to Rp E > 0.8. The numbers label the 
initial conditions of figure 2. 

FIG. 2. Values of R% E for geodesies starting from i.e. 
"topologically equivalent" to those indicated on figure 1 for 
£7=1/8, at energies from E = 0.095 up to E = 0.166. 
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